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ABSTRACT 

Real  time  simulation  of  the  time-dependent  behavior  of  low  power 
thermal  reactors  is  amenable  to  analog  computer  techniques.  This  paper 
presents  the  basic  physical  concepts  relative  to  thermal  reactor  operation, 
a  derivation  of  the  governing  transient  equations  noting  simplifying 
assumptions,  and  the  development  of  the  analogous  electrical  circuit  to 
ultimately  provide  solutions  to  these  equations.  In  addition,  many  special 
problems  involved  in  the  simulation  of  a  specific  reactor  are  discussed 
and  analyzed  in  detail.  As  a  result,  the  successful  design  and  con- 
struction of  a  simulator  whose  time-dependent  behavior  closely  approxi- 
mates that  of  an  actual  reactor  is  demonstrated.  The  reactor  simulated 
is  the  Aero  jet -General  Nucleonics  low  power  thermal  reactor  (AGN-201) 
currently  installed  at  the  United  States  Naval  Postgraduate  School, 
Monterey,  California.  Flux-time  recorder  plots  for  various  sequences  of 
operation  of  reactor  and  simulator  are  included  to  indicate  the  degree  of 
simulation  attained. 
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A  LOW  POWER  REACTOR  SIMULATOR 

1.   Introduction. 

Detailed  analysis  of  the  transient  behavior  of  nuclear  reactors 
has  been  greatly  facilitated  by  the  use  of  electrical  or  electronic 
simulators.  A  reactor  simulator  is  essentially  an  analog  computer  with 
circuitry  specifically  designed  to  solve,  on  a  true  time  base,  the 
rather  complex  differential  equations  governing  time-dependent  reactor 
operation. 

As  will  be  evident  from  the  development  of  the  governing  equations , 
complete  solution  by  direct  analysis  would  be  extremely  laborious,  if 
not  impossible,  for  other  than  a  few  limited  cases  where  modifying 
assumptions  must  be  made  before  a  useful  solution  can  be  obtained. 
Moreover,  direct  mathematical  analysis  could  not  be  expected  to  instill 
in  one  any  real  feeling  for  reactor  kinetics  under  typical  operating 
conditions.  Because  circuit  values  can  be  measured  quite  precisely, 
accurate  and  continuous  solutions  to  the  transient  equations  for  almost 
any  imposed  operating  conditions  may  be  readily  obtained  from  the  sim- 
ulator. In  addition,  by  employing  a  simple  simulator  control  system, 
a  strong  feeling  for  reactor  responses  to  any  of  a  multitude  of  typical 
operating  situations  may  be  achieved.  This  particular  feature  of  the 
simulator  might  permit  reactor  operating  personnel  or  students  to  receive 
the  equivalent  of  a  great  deal  of  reactor  control  experience  and 
operational  training  when  access  to  an  actual  reactor  is  not  feasible. 

Because  analog  computers  are  well-tried  and  accepted  engineering 
tools  for  use  by  those  involved  in  the  study  and  design  of  dynamic  systems 
of  all  types,  reactor  simulators  find  application  in  preliminary  reactor 
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design  and  especially  in  the  evaluation  and  perfection  of  associated 
systems  which  might  be  ultimately  employed  to  control  an  actual  reactor. 
In  other  words ,  completed  control  components  and  systems  may  be  connected 
to  the  simulator  by  means  of  electronic  transducers  thereby  permitting 
operational  checks  to  be  made  on  their  specific  design  functions  prior 
to  actual  installation. 

Depending  on  the  complexity  of  the  system  to  be  studied  and  the 
degree  of  simulation  desired,  simulators  may  range  in  size  from  small 
units  such  as  the  one  discussed  herein  to  extremely  large,  complex,  and 
expensive  systems  involving  the  use  of  great  numbers  of  amplifiers ,  function 
generators,  servo  and  electronic  multipliers,  electro-mechanical  devices, 
and  countless  precision  components.  In  the  case  at  hand  of  a  particular 
low  power  thermal  reactor,  utilization  of  the  basic  simulator  circuit 
together  with  a  relatively  simple  control  network  simulates  actual  reactor 
operation  to  a  surprisingly  good  degree  of  accuracy. 

The  discussion  of  basic  concepts,  the  derivation  of  governing 
equations  and  the  build-up  of  the  electrical  analog  represent  a  compilation 
of  treatments  given  on  this  subject  by  the  first  seven  authors  listed 
in  the  bibliography.  The  simulator  itself  is  essentially  a  Pagels  4 5 1 
simulator,  modified  to  duplicate  more  closely  the  responses  of  the 
individual  reactor  under  study. 

It  must  be  remembered  that  the  degree  of  simulation  attained  can  be 
no  better  than  the  degree  to  which  the  differential  equations  utilized 
define  the  system  under  study.  The  transient  equations  employed  in 
thermal  reactor  simulation  contain  many  simplifying  assumptions  and 

Numbers  in  brackets  refer  to  Bibliography  in  Section  8. 


approximations ;  and  as  a  consequence,  precise  duplication  of  actual 
reactor  behavior  may  not  be  expected.  Despite  these  assumptions,  however, 
it  will  be  seen  that  simulator  solutions  to  the  equations  employed  do 
closely  approximate  observed  reactor  behavior. 

It  is  the  purpose  of  this  thesis,  therefore,  to  present  a  rather 
complete  study  in  the  design  of  a  reactor  simulator  with  the  ultimate 
goal  being  to  construct  such  a  device  and  then  compare  time  dependent 
behavior  of  the  actual  with  that  of  the  analogous  system. 
2.   Basic  Thermal  Reactor  Concepts. 

In  order  that  the  discussion  of  the  simulation  of  a  low  power 
thermal  reactor  be  complete,  it  is  important  that  basic  concepts  be 
briefly  reviewed  prior  to  the  derivation  of  the  governing  differential 
equations . 

When  the  nucleus  of  a  fuel  atom  fissions,  it  splits  into  two  primary 
fission  fragments  of  more  or  less  equal  mass  accompanied  by  the  prompt 
emission  of  from  two  to  three  neutrons  of  high  energy.  On  the  average, 
about  2.5  neutrons  are  liberated  per  thermal  neutron  absorbed  in  fission 

c 

of  a  U23  nucleus.  The  majority  of  these  have  energies  of  about  one  or 
two  Mbv,  Most  of  the  fission  fragments  thus  formed  are  radioactive  and 
some  are  themselves  high  energy  neutron  emitters .  The  decay  constants 

for  six  delayed  neutron  precursors  have  been  measured  although  the  precise 
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origin  of  only  three  of  the  six  is  currently  known,  namely,  BR  ,  BR  ,  and 

i   .  The  average  energy  of  the  delayed  neutrons  is  about  ,5  Msv.  The 
properties  of  the  delayed  neutrons  play  a  vital  role  in  time -de pendent 
reactor  behavior  and  permit  safe  control  of  nuclear  reactor  fission 
processes  which  might  otherwise  be  extremely  difficult. 

In  addition  to  having  small  thermal  neutron  absorption  cross-sections, 


good  moderator  materials  are  those  made  up  of  elements  of  low  mass  number, 
where  the  energy  transfer  from  neutron  to  moderator  nucleus  on  collision 
may  approach  100^  as  in  the  case  of  a  hydrogenous  moderating  medium. 
This  results  in  extremely  rapid  slowing  down  of  the  high  energy  neutrons 
to  so-called  thermal  energies  (about  .025  ev) •  At  this  point,  the  neutrons 
are  essentially  in  thermal  equilibrium  with  the  moderator  atoms  (or  mole- 
cules) and  hence  are  referred  to  as  thermal  neutrons.  Despite  the 
extremely  short  slowing  down  time,  some  of  the  neutrons  will  either 
escape  from  the  reactor  or  be  absorbed  by  the  core  material  before 
thermal ization  so  as  to  be  unavailable  for  subsequent  fissioning  or  neutron 
regeneration , 

Once  the  fast  neutrons  have  been  reduced  to  thermal  energies,  they 
may  leak  or  diffuse  out  of  the  reactor  or  be  absorbed  by  moderator,  fuel, 
or  poison  nuclei! .  A  large  fraction,  about  &5%9   of  the  thermal  neutrons 
absorbed  by  the  U23  fuel  will  cause  fission;  and  a  new  neutron  lifetime 
cycle  begins  again.  The  next  job,  therefore,  will  be  to  express  in  equation 
form,  the  lives  of  neutrons  as  described  above. 
3.   Derivation  of  the  Transient  Equations. 

Based  on  the  conservation  or  balance  of  neutrons  principle,  applicable 
to  any  reactor,  a  general  expression  for  the  rate  of  change  of  neutron 
population  density  (n)  with  time  can  be  written  as  follows; 

Product ion- Absorption-Leakage  =  — cr-  (l) 

The  rate  of  production  of  neutrons  may  te  obtained  from  empirical  con- 
siderations. Likewise,  the  rate  of  absorption  depends  on  various  absorp- 
tion cross-sections  which  must  be  determined  experimentally.  For  determin- 
ing leakage,  resort  is  made  to  classical  mechanics  and  the  familiar  laws 
of  diffusion  theory  which,  although  somewhat  more  complicated,  can  be 


reduced  to  simplified  expressions  when  certain  prudent  assumptions  are 
introduced .  Let  us  now  form  the  basis  on  which  the  various  simplifying 
assumptions  can  be  made . 

The  thermal  reactor  to  be  considered  here  is  of  the  homogeneous  type 
with  reflector.  The  fuel  is  a  homogeneous  mixture  of  U23  in  the  form  of  20% 
enriched  DD2  with  polyethylene,  (CS^  ,  as  the  moderating  material.   A.  thin 
high  density  graphite  reflector  is  employed.  However,  in  the  analysis 
which  follows,  it  will  be  assumed  that  the  reactor  is  homogeneous  but  un- 
reflected.  The  degree  to  which  this  assumption  is  valid  can  best  be 
appraised  by  a  comparison  between  final  simulator  solution  and  actual 
time-dependent  reactor  test  runs.  As  will  be  seen  later,  the  assumption 
was  quite  reasonable  in  this  case . 

Before  continuing  with  additional  assumptions,  it  will  be  helpful  to 
rewrite  equation  (1)  as  the  familiar  diffusion  equation  for  thermal 
neutrons  in  non-equilibrium,  basic  to  reactor  theory,  namely 

S  -  Z.*  -   (-DV*0)  =  fg  (2) 


a 
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where  S  is  the  production  term  in  neutrons  per  cm3  per  second ,  Z   is  the 
macroscopic  absorption  cross-section  for  monoenergetic  (or  thermal)  neutrons 
in  cm  ,0  is  the  thermal  neutron  flux  in  neutrons  per  cm2  per  sec ,  and 
(-Dv*2</0  is  the  neutron  leakage  per  cm3  per  sec.  By  the  acceptance  of 
equation  (2)  as  the  equation  representing  (l),  it  has  been  assumed  that  we  are 
dealing  only  with  monoenergetic  (or  thermal)  neutrons  and  then  only  at  dis- 
tances greater  than  two  or  three  mean  free  paths  from  strong  absorbers, 
sources,  or  boundaries.  This  also  assumes  that  the  moderating  or  diffusion 
medium  is  a  poor  absorber  relative  to  its  ability  to  scatter  and,  hence,  slow 
down  the  high  energy  neutrons.  The  partial  derivative  of  n  is  utilized  because 
n  is  a  function  of  both  position  and  time.  However,  it  will  be  presumed 


that  the  reactor  is  at  or  near  critical  and  the  shape  of  the  flux 
distribution  remains  constant  even  though  it  changes  in  amplitude  with 
time.  Under  the  latter  provision,  (-DV20)  could  be  represented  by  some 
constant  multiplied  by  the  flux  0.  J-feking  use  of  this,  therefore,  let 

(-DVa0)  =2  0>  where  the  symbol  2  will  be  used  throughout  this 

1 
development  to  indicate  the  macroscopic  eross-section  for  or  the  probabil- 
ity of  occurrence  of  the  particular  process  indicated  by  the  subscript, 
i.e.,  in  this  case  Z   ■*  to  mean  probability  of  loss  by  leakage  or  leakage 
coefficient.  This  essentially  eliminates  the  space  variable  and  permits 
the  use  of  the  total  derivative  of  n,  or  tt  in  lieu  of  rr   The  flux  0 
under  the  above  considerations  is  merely  the  product  of  n,  the  density 
of  thermal  neutrons  in  neutrons  per  cubic  centimeter,  and  v,  the  mean 
thermal  neutron  speed.   Or,  expressed  in  equation  form, 

<p    =  nv   neutrons  /cm2/sec 

Thermal  neutron  flux  0,  therefore,  is  a  scalar  quantity.  We  are  now  in 
a  position  to  examine  each  of  the  terms  of  equation  (2)  in  an  attempt  to 
develop  a  more  workable  and  meaningful  equation. 

The  source  term  (S)  must  be  broken  into  several  parts,  i.e.,  S  ,  S , , 
and  Sq ,  referring  respectively  to  prompt,  delayed,  and  those  uniformly 
distributed  as  a  consequence  of  an  extraneous  outside  source  (e.g.,  a 
Ra-Be  source).  Expressed  in  equation  form: 

S  =  Sp  *  Sd  ♦  So  (3) 

The  rate  of  fission  can  be  expressed  as   2-0.  If  "V  represents  the  total 
fission  yield  of  fast  neutrons  from  both  delayed  and  prompt  sources  and 
p  represents  the  total  fraction  of  these  neutrons  which  are  delayed,  then 


the  prompt  neutron  generation  rate  is  given  by 

However,  not  all  of  these  neutrons  reach  thermal  energies.  As  previously 
mentioned  some  leak  out  or  escape  and  others  may  be  captured  by  fuel  or 
other  nuclei!  during  the  moderating  or  slowing  down  process .  Probabilities 
have  been  evaluated  for  each  of  these  processes .  The  former  is  given  by 
the  fast  non-leakage  probability,   e      (where  B2  is  the  so-called 
buckling  factor  related  to  size  and  T  is  the  Fermi  Age,  a  measure  of  the 
distance  neutrons  travel  in  slowing  from  fission  to  thermal  energies),  and 
the  latter  by  p,  the  resonance  escape  probability.  For  simplicity,  these 
factors  will  be  combined  into  a  single  thermalization  probability  symbolized 
by  P+u«  Therefore,  the  prompt  thermal  source  becomes 

Sp  =  *Zg   (i-p)-/Pth  tt) 

— B2t 

where  p..   =  pe     .  Now  the  remainder  of  the  neutrons  are  the  so- 
called  delayed  neutrons  associated  with  the  particular  fission  fragments 
mentioned  in  section  2.  The  delayed  neutron  generation  rate,  therefore, 
is  given  by 

Now,  as  previously  stated,  six  delayed  neutron  precursors  have  been 
determined  so  that  the  delayed  neutron  generation  rate  of  each  of  the  six 
may  be  expressed  by 

gi  =  ^f0  Pi 
where  p.   =  fraction  of  total  number  of  neutrons  contained  in  the  i 


delayed  group.  It  should  be  noted  here  that   B  =  ^*  B  •  At  any  time, 

1=1  i# 

there  will  be  a  concentration  C.  of  the  i   group  of  delayed  neutron 
emitters  present  in  a  unit  volume  of  the  core .  Because  these  delayed 
neutrons  are  emitted  as  a  result  of  the  norma]  radioactive  decay  of  the 
precursors,  each  emitter  obeys  the  familiar  radioactive  decay  law.  In 
other  words,  the  rate  of  change  of  the  concentration  of  the  i   group  will 
be  equal  to  the  production  rate  minus  the  decay  rate ,  or 

dCi 

dT  =   gi  "  Vi 

where  X  is  the  decay  constant  for  the  i   emitter. 
Therefore , 

dC 


-  =   Z^V"-  X1C1  (5) 


dt 


Delayed  neutrons  are  therefore  supplied  to  the  reactor  at  the  same  rate  at 
which  decay  of  the  precursors  takes  place.  This  can  be  expressed  as  a 

summation  ^      \^  C . ,   Because  the  average  energies  of  the  delayed 

neutrons  are  somewhat  less  than  those  of  the  prompt,  there  will  be  less 
tendency, for  the  delayed  neutrons  to  leak  out  during  the  slowing  down 
process^  and  hence  they  will  have  larger  values  of  non-leakage  probAWLllties. 
The   difference  in  the  values  of  non-leakage  probability  for  the  two  types 
may  be  taken  into  account  by  utilizing  effective  values  of  B..  in  the 
final  equations  which  are  slightly  larger  than  the  actual  values 0  This, 
then,  permits  multiplication  of  the  delayed  neutron  term  by  the  same  value 
of  non-leakage  probability  as  that  used  for  the  prompt  neutron  term. 
Combining  this  factor  again  with  the  resonance  escape  probability  p, 
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the  rate  at  which  delayed  neutrons  are  supplied  can  now  be  expressed  as 
a  summation  of  terms ,  namely 

i=i 


sa  =  pth 


(6) 


where  \,C.  is  defined  by  equation  (5).  Finally,  Sq  comes  from  an 
essentially  constant  extraneous  source  of  neutrons  which  makes  a  contribution 
to  production  rate  in  the  form  of  a  uniformly  distributed  source.  This 
term  will  be  discussed  in  more  detail  later. 

The   absorption  term  in  equation  (2)  needs  no  additional  explanation 
except  to  say  that  Z       represents  the  macroscopic  absorption  cross- 
section  for  all  processes,  that  is  absorption  of  thermal  neutrons  by 
structure,  moderator,  fuel,  or  poison  nuclei! . 

The  leakage  term  was  discussed  earlier  and  can  be  taken  to  equal 
2  -><t>t     where  2t   is  dependent  on  the  moderator  characteristic  and 
reactor  size. 

The  component  terms  of  equation  (2)  can  now  be  assembled  as  follows; 


dn 
dt 


P  d 


So)    -  2a  <fi  -   Z1  <p 


dn     _ 
dt 


*Zf(l-  p)Vpth  +Pth 


\ci  ♦  S0 


or 


-♦(2a*2x) 

fs  =  to2t[i-%f  *;|  Vl)PJ    -*(2a*2l)  *  s„       (7) 


La  order  to  reduce  equations  (5)  and  (7)  to  more  usable  forms,  several 
terms  must  be  defined.  First,  the  term  reactivity  jO   may  be  defined  as 


(?) 


eff 


where  K  ~f  is  the  effective  multiplication  factor  for  a  given  finite  sized 
reactor.  By  effective  multiplication  factor  is  meant  the  population  ratio 
of  successive  neutron  generations.  If  1  is  next  defined  as  the  mean  lifetime 
of  thermal  neutrons  in  the  reactor  and  1/Iq  as  the  rate  at  which  neutrons , 
prompt  and  delayed,  are  thermalized  in  the  reactor,  then  K  «f  could  also 


be  expressed  as  the  ratio  of  1  /1q      56  J  , 


or 


Keff  =   1*  A*0  (9) 

This  may  be  interpreted  as  follows:  K  _«  is  the  ratio  of  the  time  spent 
by  the  thermal  neutrons  diffusing  in  the  reactor  to  the  time  spent  in  being 
thermalized.  In  other  words,  if  it  takes  longer  to  thermalize  the  second 
generation  of  thermal  neutrons  than  the  life  expectancy  of  the  first 
generation,  the  chain  can  be  expected  to  die  out.   In  this  case,  K  ff  would 
be  less  than  unity.  The  reverse  would  be  true  if  K  „f  were  greater  than 
unity.  When  K .  - «  is  precisely  unity  (or  in  the  case  of  the  presence  of 
an  extraneous  source,  slightly  less  than  unity),  just  as  many  thermal 
neutrons  are  born  as  die  and  the  chain  is  sustained.  These  three  situations 
represent  subcritical,  supercritical,  and  critical  behavior  of  the  reactor. 
10  can  be  expressed  as  the  reciprocal  of  the  t normalization  rate  J  6  t  ,  or 

*  1 

*o  =   ^ — 7 (10) 

vZfVpth 

This  term  is  also  given  the  name  of  modified  lifetime  by  Murray  ?  3 j  • 
Now  thermal  neutron  lifetime  in  an  infinite  medium  would  be  given  by  the 
average  distance  traveled  by  the  neutron  before  absorption,  called  the 


absorption  mean  free  path,   X  (  =  — = — ),  divided  by  the  mean  neutron 

a     2?a 
velocity  v.  However,  in  a  finite  reactor,  this  must  be  reduced  by  the 


10 


thermal  leakage  factor 


Thermal  Absorption 


Thermal  Absorption  ♦  Thermal  Leakage 


or 


'Z     2a  +  Si     y(2,*2i) 


(11) 


Substituting  (9),  (10),  and  (11)  into  (?),  therefore, 


/°  = 


1    - 


iX 


=  1  - 


T(2t*Z!) 
v2f"^Pth 


(12) 


.«  .* 


Now,  utilizing  1  ,  10  and  /O   as  defined  above,  equations  (5)  and  (7)  mar 
be  simplified  as  follows: 


dn 
dt 


Lzf  [i-P}V>  gvi^th 


-  *  (2a  *Zl)    ♦  s, 


ff     =       |-f  n   (1  -  p)   ♦  p 


tb    -^ 


_    n 


♦     S, 


■S-  (  1  -  p  -         ■      )   +  p        J^L  X.C.     +    s, 

LS  Keff  th     i=l     i  1  ( 


or  it  =  ^t-8-  n  ♦  Pth 


i=i 


x1ci    ♦    s0 


(13) 


dC, 


Likewise,  from   (5),       ^     =     ( Zf  0^"/"-    J^     > 


dC, 


or       \ci     = 


^  Pth 


n     - 


dt 


(H) 


Equations  (13)  and  (14-)  are  the  desired  transient  equations. 
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Finally,  a  single  governing  equation  can  be  written  by  a  simple 
substitution  of  (14)  into  (13)  as  follows: 


dt    '  ,*        Fth  *£-.  x  ,*  dt 


1?  i=i    lo  P 


th 


6  dC 

=  /°   n  -  — *-  n  +    |    n  -  p    ^>  t~  + 

/  *        ,*        -,*        th  f~  dt 
lo        ^o        -^           i-1 


So 


■  jj~  n  *  pth  2   dT  * 


S0 


-o  i=l 


-<*  ■/•■-.\i4.<» 


It  will  be  further  assumed  that  the  so-called  one-group  theory  is 
sufficiently  accurate  for  accomplishing  the  goal  at  hand.  In  one  group 
theory  it  is  supposed  that  all  production,  diffusion,  and  absorption  of 
neutrons  occur  at  a  single  energy,  namely,  the  thermal  energy.  Although 
this  is  a  considerable  approximation,  the  final  simulator  circuit  is 
thereby  simplified  and  results  are  still  quite  satisfactory.  Since  this 
theory  implies  that  all  fission  neutrons  have  thermal  energies,  the 

resonance  escape  probability,  p,  and  the  non-leakage  probability  during 

-B2T 

the  slowing  down  process,  e      ,  do  not  affect  the  multiplication  factor 

and,  hence,  become  unity  in  equations  (13)  thru  (15).  The  two'  terms  are 
carried  through  in  the  development  of  the  electrical  analog  but  no 
attempt  is  made  to  employ  thera  as  variables  in  the  final  circuit. 
4,   The  Electrical  Analog. 

Equation  (15)  is  now  in  a  convenient  form  for  a  computer  solution. 
The  network  given  in  Figure  1  will  be  shown  to  describe  this  differential 
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equation.  Referring  to  Figure  1,  and  using  Kirchoff's  second  law,  it 
is  evident  that 

6 

h  =  h  -  IE    h  +   x9 

i=l 

Note  the  similarity  in  form  to  (15). 

(1*  £s   \ 
The  Rate  Term  v  °  dt  ; 


t         32 
ll  c7  dt 


and  letting  E  in  volts  be  proportional  to  n  in  neutrons  per  cm3  per 
sec  by  the  constant  U  ,  i.e.,  E  =  o  n,  then 

3   c7  1   .*   dn 


■HH 


I7  =   "^M  lo   |?  (16) 


The  Reactivity  Term  COTl) 

Let  A  be  the  proportional  distance  of  the  potentiometer  tap  (at  P) 
from  the  center  of  the  potentiometer.  Referring  to  the  more  detailed 
development  given  in  Appendix  I, 


T  AE 

*8 


R  [l  *  B/P     ♦  PAR  (1  -  A2F] 


The  range  of  values  of  A  are  zero  at  the  center,  +  1  at  the  ♦  E  end, 
and  -1  at  the  -E  end.  Now,  it  is  apparent  that  L>  will  have  its 

maximum  positive  value  at  A  =  ♦  1.  Therefore,  the  range  of  reactivity 

variation  of  the  simulator  may  be  set  and  the  proper  scaling  factor 

determined  by  setting 
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I, 


8 


/^nax  c  max 


or        /O    = 


/        i  ♦      pa 


4fl(R  ♦  P)     (1  '  A2) 


Since  L*      occurs  at  A  =  1,  then 
o  max 


E 
h   max      R(l  +  R/P) 


and  therefore, 


Aax  <17> 


X8       R(l  +  R/P)     7^ 


max 


\gain  representing  E  by  Jn,  this  reduces  to 


I*   =   7  p/i  :  pa>s  >n — (/?*  (1?) 


L8   "   /  "RlT 


***Y~}>f n 


More  careful  consideration  of  (17)  above  reveals  that  /O   is  a  non-linear 
function  of  A  which  is  undesirable  for  our  purposes .  This  is  due  to 
the  potentiometer  loading  having  its  greatest  effect  around  the  mid-position 
of  the  potentiometer.  Linearity  can  be  very  closely  approximated ,  how- 
ever, by  making  R  sufficiently  large  in  relation  to  P  that  the 
denominator  of  (17)  becomes  essentially  unity,  leaving 

f  ~     max 
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If,  for  example,  R=5P»  the  linearity  remains  below  one  percent  for  any 

value  of  A  from  (-)l  to  (-01. 

,  *      ,6   dC. 

(1°  pth^dr} 

The  Delayed  Precursor  Term i=l 

6 
The   ,^!1  term  of  (15)  is  simulated  by  the  six  R.c.  delay  networks 

The  basic  voltage  equation  describing  the  voltage  drop  across  any  of  the 
i    type  networks  is 


Qi  dQi 

E  =  —r-     +   R.   tt1 
c„        i   dt 


where  Q.   is  the  capacitor  charge  in  coulombs  and  c.   is  the  capacitance 
of  the  i    group.  The  electrical  charge  is  analogous  to  the  "storing" 
of  delayed  neutrons  in  the  fission  fragments  before  they  are  emitted.  The 
corresponding  delayed  neutron  equation  represented  by  these  networks  is 
(14.) ,  which  upon  rearranging  becomes 

*•  Pth  xici  ^  Pth     dci 

n  —        "        ♦   — — — 

Pi  Pi     dt 

Again  letting  E  =0  n  and  matching  terms  with  the  basic  equation  above, 
3   Xo  pth  1 

Ql  =     hh      c±  and  Ri Cl  =  ^T 


A    ]    lo  Pth    dT 


Now«     h   '  it-   =   ( -jrr  >   ^  Pth  ^  ; 


dC. 
and  by  substitution  for  gr~  from  (14) , 
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h  =  <  -$q- >  h n  -  ( -tjsr }  x*0  Pth  xici 


But,       ^     I.      =     I|    ♦  Ia  ♦  I3   .    .    .    . 
=1       X 


and 


<  pth  £  =  j^fe-}  ^ n "  [l^-]15^  f?  Vi  (19) 


?iRl 
Therefore,        p.R.      =     fcRi    =     p2R2     =     P3R3 

The  Source  Term   (S.) 

If    Rg>>  P9,     then 


9    '^        Rg 


where  Ag  is  the  fraction  off  ground  potential  and  V  is  a  constant 

So. 

(max) 


voltage  source.  Now,  letting  AQ  =   Q  ^ then 


.* 


Iq  =   j    l —  (   lo  S0  (20) 

*9  S°  (max)10 


Finally,  summing  the  currents  as  given  by  equations  (16),  (1?),  (19) 
and  (20), 

6 
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^      7 1  ,*     dn 


■w ♦  ^Hajr  -  y&r$ 


0  Pth  f=r    «~ 


{ 


V^max)  ^ 


1*  so  (21) 


In  order  for  the  current  output  of  each  term  to  be  of  the  correct 
magnitude,  the  scaling  factors  which  are  enclosed  by  brackets  in 
equation  (21)  must  be  identically  equal  to  each  other;  that  is, 

c-$"i =  [*<i  ♦&>/>■«} =  Hfc? =  [^(L)1-]      (22) 

If  this  can  be  done,  then  these  coefficients  drop  out  of  equation  (21) 

leaving  equation  (15)  which  is  to  be  duplicated  electrically.  Values  to 

insert  into  (22)  will  be  determined  by  constants  in  the  actual  system  to 

be  simulated  as  well  as  arbitrary  values  selected  in  accordance  with  standard 

computer  operating  techniques. 

5.   Simulation  of  a  Specific  Low  Power  Thermal  Reactor  — AGN-201 

General  thermal  reactor  simulation,  therefore,  can  be  readily  duplicated 
by  the  circuit  given  in  Figure  1,  This  circuit  when  set  up  on  a  conven- 
tional analog  computer,  using  completely  arbitrary  circuit  parameters,  will 
yield  typical  time-flux  curves.  The  problem  at  hand,  however,  is  to 
simulate  a  specific reactor  wherein  lie  many  additional  considerations  which 
first  must  be  resolved  before  a  final  circuit  can  be  drawn  up.  These 
include  detailed  investigations  to. 

(1)  determine  actual  range  of  reactivity  variation,  compute  rates 
of  reactivity  change  for  control  rod  movement,  and  provide  suitable  means 
for  duplicating  control  rod  action; 

18 


(2)  determine  a  means  of  introducing  a  maximum  reactivity  control 
and  temperature  coefficient  of  reactivity; 

(3)  solve  electrical  problems  resulting  from  real  time  simulation  over 
wide  operating  ranges  of  neutron  flux  from  source  level  to  full  power j 

(4)  and  to  determine  an  appropriate  value  of  neutron  source  voltage 
corresponding  to  the  extraneous  source  production  rate. 

In  the  analysis  of  these  problems,  attention  shall  be  restricted  to 
the  simulation  of  the  low  power  thermal  reactor  currently  installed  at  the 
United  States  Naval  Postgraduate  School.  This  reactor  will  hereinafter  be 
referred  to  as  the  AGN-201. 

Before  solving  the  problems,  it  is  necessary  to  mention  a  few  essential 
features  of  the  AGN-201  which  will  form  the  basis  on  which  ultimate 
simulator  construction  may  proceed.  The  general  characteristics  of  the  201 
reactor  are  listed  in  Table  1  below. 

TABLE  1.  GENERAL  CHARACTERISTICS  AGN-201 


Power  (licensed) 

Maximum  Thermal  Flux  @  100  mw 

Critical  Mass 

Core  Size  (cylinder) 

Maderator  (polyethylene ) 

Reflector  (graphite) 

Shield  (lead) 

(water) 
Over-all  diameter 
Over-all  height 

Temperature  coefficient  of  reactivity 
(based  on  water  shield  teap.  data) 


100  miliwatts 

4,5  x  10  cm"2  sec"" 

665  gm  U2*5 

10"  dia.  x  9£»  high 

11  kgm 

20  cm  thick  (1.70  gm/cm3) 

10  cm  thick 

55  cm  thick  (1000  gallons) 

6J  ft 

9  ft 

-2.75  x  10"4  />C 
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Figure  2  shows  the  A.GN-201  core  tank  containing  the  uranium-polyethylene 
core  and  portions  of  the  graphite  reflector.  Figure  3  is  a  cross-sectional 
view  of  the  entire  reactor  unit  showing  the  core,  reflector,  lead  and 
water  shields,  glory  hole,  and  control  rod  assemblies.  The  critical  mass 
consists  of  about  6A5  grams  of  U23  in  the  form  of  20$g-  enriched  U02.  It 
is  contained  in  the  form  of  a  powder  dispersed  in  a  solid  polyethylene 
moderator.  The  core  is  made  in  nine  sections  of  various  thicknesses 
(Figure  2).  Each  of  the  four  bottom  discs  has  three  5-cm  holes  and  one 
2.5-cm  hole.  The  201  has  two  safety  rods  and  two  control  rods.  The 
three  larger  holes  are  for  two  safety  rods  and  one  coarse  control  rod  while 
the  smaller  hole  is  for  a  fine  control  rod.  The   two  safety  rods  and  the 
coarse  control  rod  are  identical  and  each  contain  about  14-  grams  of  U23  . 
The  fine  rod  containaabout  3.5  grams  of  U23  .  All  rods  operate  so  that 
reactivity  is  increased  as  the  rod  is  inserted.  The  amount  of  reactivity 
each  rod  controls  increases  as  the  amount  of  active  material  contained 
increases.  With  the  above  amounts  of  U23  in  the  rods,  the  safety  and 
coarse  control  rods  control  about  1.3$  reactivity  while  the  fine  control 
rod  controls  about  0 .%   reactivity.  The  rods  are  inserted  by  a  reversible 
d-c  motor  drive  mechanism  requiring  about  one  minute  to  travel  the  full 
25-cm  stroke.  Referring  to  Table  1,  the  thermal  flux  at  licensed  full 
power  is  given  as  U, 5  x  10  neutrons  /cm2/sec#  This  is  well  below  that 

at  which  Xenon  poisoning  becomes  a  consideration,  namely,  at  a  flux  in 

12 
excess  of  about  10   neutrons  /cm2/sec.  Finally,  because  of  the  low 

power  developed,  heating  of  the  core  is  so  slight  that  it  is  essentially 

ambient  temperature  which  determines  core  temperature  and,  hence,  the 

maximum  reactivity  available .  The  negative  temperature  coefficient  listed 

in  Table  1,  therefore,  is  based  on  the  change  in  temperature  of  the  1000 
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gallon  water  shield  surrounding  the  core  tank  i  10  (  ,   With  these  features 
in  mind,  let  us  proceed  with  the  investigations  into  the  special  problems 
mentioned . 

First,  a  control  rod  calibration  must  be  made,  the  result  of  which  is 
shown  graphically  in  Figure  A,     The  term  relative  reactivity  as  used  here 
implies  an  arbitrary  zero  reactivity  reference  when  the  rods  are  fully 
retracted.  As  will  be  seen  later,  the  major  portion  of  the  reactivity 
contribution  or  "worth"  of  the  rods  is  employed  in  reducing  the  amount 
by  which  the  reactor  is  subcritical.  Only  a  small  fraction  of  this 
relative  reactivity  is  actually  positive  reactivity  as  defined  by  (8). 
This  curve  represents  a  superposition  of  data  taken  on  both  fine  and 
coarse  control  rods  which  actually  exist  on  the  reactor,  thus  providing 
a  single  curve  to  represent  the  movement  of  both  rods  simultaneously. 
In  view  of  the  fact  that  the  action  of  the  two  control  rods  on  the  reactor 
can  be  duplicated  more  simply  and  more  economically  on  the  simulator  by 
a  single  controlling  device,  simulation  of  operation  of  each  rod  individ- 
ually does  not  seem  warranted.  Now,  under  normal  conditions,  the  reactor 
will  be  operating  in  the  straight  line  portion  of  the  curve  when  near- 
critical  conditions  exist.  As  indicated  in  Figure  4.,  the  slope  of  this 
straight  line  is  ?37  \isO   Lox   10  )  per  cm  of  rod  travel.  From  weight 
data  on  the  fuel  discs  which  make  up  the  two  safety  rods  and  the  fine 
and  coarse  control  rods,  it  is  found  that  the  relative  reactivity  con- 
tribution of  each  safety  rod  is  essentially  equal  to  that  of  the  coarse 
control  rod  or  about  11,275  p/O  each.  Therefore,  the  total  contribution 
in  terms  of  relative  reactivity  of  all  rods  is  found  to  be  36,725  p/o    • 
The  reactivity  potentiometer  can  now  be  designed  to  match  the  control  and 
safety  rod  curves.  Figure  5  shows  a  plot  of  percent  of  total  resistance 

t 
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of  the  reactivity  potentiometer  versus  percent  of  total  rotation  (about 
350  )  where  all  rods  are  being  considered.     Notice  that  the   curve  is 
shown  up  to  the  point  at  which  the  safety  rods  are  completely  inserted 
as  a  dashed  line.      Because  the  safety  rods  are  not  utilized  for  control 
purposes,   their  individual  calibration  curves  are  of  no  concern  as  long 
as  their  ultimate  contribution  to  total  relative  reactivity  is  known, 
namely,   22,550     y>p      or  about  61. i$  of  the  total  available  for  addition 
by  rod  insertion.     From  this  point  on,   however,   it  would  be  desirable  to 
follow  the  remainder  of  the  curve  which  is  essentially  Figure  4  super- 
imposed but  plotted  using  non-dimensional  relationships.     To  duplicate  the 
entire  curve  precisely  would  require  a  non-linear  potentiometer  which  could 
be  designed  but  at   considerable  expense.     Instead,   in  view  of  the  fact  that 
the  straight  line  portion  of  the  combined  control  and  calibration  curve 
is  of  greatest  interest,   a  center-tap  potentiometer  permits  precise  duplica- 
tion in  this  region.     At  5C#  rotation,   extrapolation  of  the  straight  line 
portion  back  to  this  point  reveals  that  25^  total  resistance  should  result 
at  this  point.     For  the  20,000  ohm  potentiometer  selected,   this  requires  a 
5000  ohm  resistor  in  parallel  with  the  first  half  of  the  potentiometer.     The 
arrangement  is  shown  schematically  in  Figure  5;   and  the  actual  plot  of 
percent  resistance  versus  percent  rotation  resulting  is  indicated  by 
circles  in  the  figure.     It  is  noted  that  61.4^  of  the  total  resistance 
will  now  occur  at  about  75/£  of  total  potentiometer  rotation.     If  this  point 
is  denoted  for  the  simulator  as  having  a  control  rod  position  reading  of 
zero,   then  a  discrepancy  will  exist  between  the  control  rod- re activity 
relationships  of  the  actual  reactor  and  that  of  the  simulator.     As  seen 
from  Figure  5,   the  simulator  rods  increase  reactivity  linearly  over  their 
entire  travel  while  actual  rod  insertion  increases  reactivity  linearly 
only  over  the  last  10  cm  of  travel.     This  situation  must  be  accepted  when 
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the  actual  reactivity  curve  cannot  be  duplicated  over  its  entire  range. 
The  rate  of  potentiometer  rotation  is  easily  obtainable  by  timing  the 
actual  control  rod  insertion  rate.  Both  fine  and  coarse  control  rods  have 
been  set  to  move  from  the  fully  retracted  to  the  fully  inserted  positions 
(0  to  25  cm.)  in  about  one  minute  at  the  fast  speed  setting.  Therefore, 
the  last  10  cm  of  actual  rod  travel  will  be  completed  in  0.40  min.  From 
Figure  5  the  last  10  cm  of  actual  rod  travel  commences  at  85%  of  total 
potentiometer  rotation  or,  alternately,  over  the  last  0.15  revolution. 
Hence,  the  reactivity  potentiometer  drive  rate  should  be  0.15  revs/0,40 
min  =  .37  rpm,  at  the  proper  rate.  Over  this  straight  line  portion, 
reactivity  change  at  fast  speed  is  found  to  be 
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^jo   /cm  x  7q   cm/sec  =  349  v>£>   /sec 


To  solve  the  second  problem,  a  means  of  establishing  a  variable  zero 
or  critical  reference  level  must  be  provided  for  the  reactivity  potentio- 
meter. The  reactor  functions  between  certain  definite  reactivity  limits 
for  any  specified  core  temperature.  In  other  words,  some  scheme  is  re- 
quired to  duplicate  the  normal  operating  ranges  to  be  expected  from 
maximum  subcritical  conditions  to  maximum  supercritical  conditions.  For 
example,  at  23°Cf  the  maximum  positive  reactivity  attained  with  all  rods 
fully  inserted  is  only  about  .0020.  At  the  same  time,  with  all  rods  fully 
retracted,  the  reactor  has  a  negative  reactivity  of  about  .035.  This 
one-sided  variation  is  representative  of  the  type  of  reactivity  behavior 
encountered  under  typical  operating  conditions.  It  should  be  noted  here 
that  the  above  values  were  determined  with  the  so-called  "glory  hole"  empty. 
Ike  glory  hole  (Figure  2)  is  a  small  access  tube  about  one  inch  in  diameter 
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which  passes  through  the  core  center,  and  into  which  additional  fuel, 
moderator,  reflector  material,  samples  and  the  like  can  be  inserted  to 
yield  changes  in  reactivity.  As  an  example,  by  insertion  of  the  proper 
amounts  of  moderator,  reflector,  and  lead  shielding  in  the  glory-hole, 
the  core  is  made  essentially  homogeneous;  and  an  additional  reactivity 
contribution  of  about  2600  y>p      can  be  realized.  Moreover,  an  average 
linear  temperature  coefficient  of  reactivity,  based  on  water  shield  tem- 
perature, of  about  -2o75  x  1CH*  per  °C  exists  for  the  AGN-201.  This  rather 
significant  effect  on  reactivity  must  also  be  taken  into  consideration. 
Figure  6  represents  the  scheme  devised  to  yield  the  proper  operating  range 
of  reactivity  from  a  relatively  large  negative  value  to  a  small  positive 
value  and,  in  addition,  to  permit  control  over  maximum  positive  reactivity 
available.  A  permanently  imposed  off balance  in  the  values  of  R,   and  R2 
will  electrically  offset  the  mechanical  mid-position  as  indicated.  The 
values  given  in  the  figure  are  those  based  on  the  reactivity  conditions 
at  23 °C  and  employing  a  homogeneous  core0  The  critical  or  zero  position 
is  now  located  at  point  C.  In  addition,  it  is  noted  that  C  can  be  further 
controlled  by  movement  of  the  maximum  reactivity  control  and  temperature 
coefficient  (MRCTC)  potentiometer.  If,  for  example,  the  temperature  in- 
creases (or  alternately,  if  a  sample  with  a  large  thermal  neutron  absorption 
cross-section  is  placed  in  the  glory  hole),  the  MRCTC  potentiometer  is  moved 
as  indicated  in  the  figure.  This  increases  the  resistance  in  loop  2  and 
decreases  the  resistance  in  loop  1,  causing  an  unbalance  in  the  feedback 
circuits,  thereby  producing  a  decrease  in  reactivity  corresponding  to  the 
amount  of  resistance  change  applied.  If,  therefore,  it  is  desired  to  hold 
the  reactivity  at  zero,  the  reactivity  potentiometer  must  be  moved  in  such 
a  direction  as  to  balance  the  resistances  in  the  two  loops.  In  this  case, 
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it  would  require  movement  of  the  arms  as  indicated.  It  can  readily  be  seen 
that  if  the  temperature  coefficient  is  offset  far  enough,  the  maximum 
positive  reactivity  available  can  be  reduced  until  the  reactor  may  no  longer 
be  brought    critical  even  though  all  rods  are  inserted.  This  corresponds 
to  the  point  marked  B  in  the  figure.  The  original  electrical  center  position 
C,  therefore,  effectively  walks  up  and  down  the  reactivity  potentiometer 
depending  on  the  setting  of  the  MRCTC  potentiometer.  Of  course,  maximum 
positive  reactivity  may  likewise  be  increased  to  prompt  critical  conditions 
and  above  if  desired.  One  of  the  solutions  of  the  transient  equations  for  a 
step  change  in  reactivity  is  the  so-called  Inhour  equation.  The  plot  of  this 


equation  using  the  value  of  neutron  lifetime  applicable  to  the  AGN-201  is 
given  in  Figure  7.  Calibration  of  the  MRCTC  potentiometer  can  be  made  by 


measuring  the  stable  period  (T)  for  any  particular  potentiometer  setting 
when  all  rods  are  fully  inserted.  The  corresponding  reactivity  may  then  be 
picked  off  an  appropriate  plot  of  the  Inhour  equation  using  the  neutron  life- 
time applicable  to  the  reactor  being  simulated.  Such  a  calibration  curve  is 
given  in  Figure  8  for  the  particular  MRCTC  potentiometer  employed  here. 
Temperatures  shown  in  the  figure  are  based  on  experimental  data  employing  a 
homogeneous  core  and  using  the  negative  temperature  coefficient  of  reactivity 
from  Table  1.  In  summary,  therefore,  the  MRCTC  potentiometer  functions  as 
an  electrical  stop  on  the  upper  limit  of  the  reactivity  potentiometer  which 
provides  a  scheme  to  preset  any  desired  maximum  reactivity  (or  alternately, 
any  desired  core  temperature).  The  effect  on  the  electrical  circuit  can  be 
taken  into  account  quite  simply  as  is  indicated  in  Appendix  II. 

The  third  problem  can  be  solved  by  employing  stablized  operational 
amplifiers  with  virtually  zero  drift  to  give  reasonably  accurate  and  repro- 
ducible solutions  at  the  extremely  low  starting  or  source  voltages  encountered. 
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The  latter  are  of  the  order  of  magnitude  of  a  few  millivolts  if  scaling  is 
such  as  to  have  about  four  decades  of  multiplication  from  source  level  to 
full  power.  The  maximum  output  of  the  amplifiers  is  100  volts.  A  variable 
biasing  arrangement  is  also  desirable  to  ensure  zero  offset  at  zero  voltage 
input  to  the  amplifier  gird.  Any  one  of  several  suitable  methods  recom- 
mended by  the  manufacturer  of  the  amplifiers  will,  satisfy  this  requirement. 
Causing  additional  problems  at  low  operating  voltages  are  the  condensers  in 
the  circuit.  It  is  found  that  an  erroneous  flux-time  curve  results  in  the 
low  voltage  region  upon  the  decay  of  the  neutron  flux  when  the  reactor 
simulator  is  scrammed.  Analysis  of  oil,  paper,  mica,  electrolytic,  and 
other  types  of  commonly  used  capacitors  reveals  a  marked  dielectric  ab- 
sorption effect  by  which  these  types  hold  their  charge  for  a  relatively 
long  period  of  time  even  after  direct  shorting  across  the  terminals.  This 
problem  is  virtually  eliminated  by  the  use  of  polystyrene  or  mylar  condensers 
which  exhibit  characteristically  low  dielectric  absorption  properties. 


The  fourth  and  final  problem  is  significant  only  in  that  it  is  con- 
sidered desirable  to  simulate  reactor  responses  over  the  entire  operating 
range.  This  demands  computer  solutions  to  the  controlling  equations  at   - 
extremely  low  values  of  neutron  flux  which  corresponds  to  the  condition  in 
the  reactor  when  all  rods  are  retracted,  and  the  existing  neutron  density  is 
dependent  on  the  extraneous  neutron  source  strength.  This  problem  may  be 
approached  as  follows:  From  equation  (15),  namely, 


*  dn  .,*     -£-  dCi      * 

*•  dt    =  /°  n  -  ^  pth  2;  dF  +    *•  s< 


0 


dn       dCi 
it  is  apparent  that  under  steady— state  conditions  -rz     and   ..   terms 

drop  out  leaving 
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=  0 


/O   n  +  1q  S 

or  rearranging , 

So  =   ; (23) 

And  from  equation  (8)  and  (9)t  this  can  be  written  as 

n(l  -  K   ) 
S,  =   ^6ff  (24) 

From  manufacturer's  data  or  from  actual  experimental  flux  determinations,  a 
value  of  n  can  be  calculated  at  a  level  where  reactivity  is  known.  The 
latter  is  obtainable  from  experimentally  determined  control  rod  calibration 
data.  Employing  (23) »  this  will  establish  a  value  of  S010  at  one  particular 
control  rod  position.  From  the  known  value  of  reactivity  at  this  position, 
Keff  m8i^   k®  determined  from  (P).  Finally,  from  (24-),  SqI*  is  obtained. 
If  the  value  of  estimated  neutron  lifetime  as  given  by  the  manufacturer  is 
then  employed,  S0  is  readily  calculated.  It  is  noted  from  (22)  that 


J 


¥•  w  tf  R(1  +  ** V° 


max 


or  rearranging, 


=    J  h   (S°(mi)  *•  } 
R(l  ♦  B/P)  /^>] 


max 


Therefore,  to  obtain  the  neutron  source  voltage  corresponding  to  a  particular 
source  generation  rate,  the  product  (S0 /   x  Iq  )  need  only  be  determined 
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as  accomplished  above,  and  actually  the  two  terms  need  not  be  separated 
to  determine  the  voltage.  The   magnitude  of  neutron  lifetime  itself  is 
important  in  the  determination  of  a  proper  value  of  the  integrating 
capacitor  c7.  Prom  (22) 

C7  ~    R(l  ♦  R/P)  , 

max 


/°. 


or,  essentially,  c„  is  proportional  to  neutron  lifetime.  However,  if 
the  order  of  magnitude  of  neutron  lifetime  is  known,  this  is  sufficiently 
accurate  to  yield  excellent  results  as  long  as  maximum  reactivity  values 
are  limited  to  .005  and  under.  This  figure  of  .005  is  arrived  at  by  evalua- 
ting the  stable  period  corresponding  to  various  values  of  positive  re- 
activity when  different  neutron  lifetimes  are  considered.  The  Inhour 
Equation  may  be  employed  to  make  this  comparison  analytically f  and  it  will 
be  found  that,  at  low  reactivities,  period  is  independent  of  neutron  life- 
time for  a  specified  reactivity.  Furthermore,  in  view  of  the  fact  that  the 
stable  period  corresponding  to  a  neutron  lifetime  of  10"  sec  and  a  re- 
activity of  .005  is  about  2.5  seconds,  this  condition  on  the  actual 
reactor  cannot  be  duplicated  because  a  period  scram  is  set  to  function  at 
about  3.0  seconds.  As  a  consequence,  comparison  runs  at  higher  values  of 
reactivity  to  better  evaluate  neutron  lifetime  are  not  feasible 0 

Having  solved  the  above  four  problems,  the  remaining  circuit  para- 
meters may  be  calculated  as  indicated  in  Appendix  II  and  the  final  AGN-201 
simulator  can  be  constructed.  Figure  9  is  a  schematic  of  the  final  analog, 
and  listed  in  Table  2  below  are  the  various  values  of  circuit  components 
referred  to  in  the  figure.  Also  listed  are  the  other  factors  which  appear 
in  the  scaling  factor  equation  (22). 
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TABLE  2a.     CIRCUIT  VALUES 


p     = 

13.0  Kjl 

(actual  value  of  pot.) 

R,      = 

119     Ka 

R*    = 

129     Kii 

c,      ■ 

.0007     »if 

H, 

=    102M 

°»  = 

.0207    j*f 

*a 

=       30M 

0,      = 

.208       jif 

«3 

=    10.6M 

\  * 

.5U      pf 

\ 

*     12  .OM 

°5    = 

2.06      >tf 

*5 

=     15. AM 

c6    * 

.785      |if 

% 

=     102M 

°7    " 

.0039    »if 

Rf      =      1M 

S  s 

10M 

V      « 

51  mv 

TABIE  2b.  ADDITIONAL  SCALING  FACTOR  VALUES 

P         - 
reff 

23.0K.ft 

(use 

for  P  in  equation  22) 

R 

119KJZ 

(use 

for  R  in  equation  22) 

/°max= 

.034775  i 

23°C 

(core  homogeneous) 

Pt       c 

.00025 

Pa      * 

.00085 

Pi      * 

.002a 

»*      = 

.00213 

h    * 

.00166 

P6      " 

.00025 

J>          = 

2.22 

So       » 

S°(max) 

=     .579  neuts/cmVsec 
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1*  &  1.0  x  10"*  sec  9*   l£  near  critical 
(i.e.  Keff  &   1) 

Figure  10  is  a  photograph  of  the  completed  simulator  with  linear 
recorder  unit  attached.   Important  features  and  actual  circuitry  are 
shown  in  Figure  11  which  is  a  rear  view  of  the  simulator  console.  At 
the  left  center  of  Pigure  11  is  the  reactivity  potentiometer  and  drive 
motor  assembly |  at  the  upper  right  portion  of  the  figure  is  the  MRCTC 
potentiometer;  and  at  the  bottom  portion  of  the  photograph  is  the  entire 
simulator  circuit  which  is  constructed  to  resemble  the  schematic  of 
Figure  9,  Along  the  top  of  the  console  are  located  the  operational 
amplifiers  and  stabiliting  amplifiers. 
6.   Degree  of  Simulation  Attained 

Figures  12  through  14  represent  comparison  flux-time  plots  taken 
from  the  reactor  and  the  simulator,  utilising  the  same  linear  recorder 
so  as  to  ensure  an  identical  time  base.  The  sequence  of  events  describing 
each  of  the  comparison  runs  is  briefly  discussed  in  the  succeeding 
paragraphs . 

The  sequence  of  events  describing  Figure  12  is  given  in  Table  3. 

TABLE  3.    SEQUENCE  OF  EVENTS  FOR  FIGURE  12 

TIME  IN               REACTOR  SEQUENCE  CORRESPONDING 

MINtnSS \ SIMULATOR  SEQUENCE 

Steady  flux  level;  rods  @  15         Steady  flux  level; 

cm;  glory  hole  empty;  temp.  rods®  10.1  cm;  MRCTC 

21° C.  set  for  proper  max. 

reactivity  conditions . 

0        Inserted  rods  @  fast  speed.  Inserted  rods  @  fast 

speed. 
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0.4 


2.1 


4.25 
7.50 


8.65 


10.15 


12.05 


13.5 


Rods  up;     flux  permitted  to 
increase  on  17.5  sec.  period 


Flux  level_g  7.6  x  10 
(0.76  x  10       mia  scale  read- 
ing); rods  retracted  to 
20.4-5  cm;  flux  permitted  to 


decay  to  1.0  x  10 


Adjusted  rods  for  steady 
flux  level. 

Quickly  inserted  moderator, 
graphite,  and  lead  discs  in- 
to glory  hole  to  provide  a 
step  change  of  about  (♦) 
1300  pa   \   of  total  required 
to  make  core  essentially 
homogeneous , ) 

Retracted  rods  to  19.93 
em;  allowed  flux  to  de- 
cay to  4.0  x  106. 

Inserted  remainder  of 
poly,  graphite,  and 
lead  in  glory  hole  and 
adjusted  rods  to  steady 
flux  level  of  4.5  x  10 
?  rated  power  100  aw. 


Quickly  withdrew  all 
poly  from  glory  hole  to 
provide  a  step  change  of 
about .{-)   2600  p^  ;  al- 
lowed flux  level  to  decay 
to  about  0.63  *  106o 


Inserted  rods  all  the 
way;  permitted  flux  to 
increase  to  about  4.5 
x  10"  and  leveled  at 
about  4.0  x  106. 


Rods  up;  flux  per- 
mitted to  increase 
on  17.5  sec.  period. 


Flux  level 
,6 


7.6  x 


10v  (76  volts  scale 
reading)  rods  re- 
tracted to  22.5  cmf 
flux  permitted  to^ 
decay  to  1,0  x  10  • 

Adjusted  rods  for 
steady  flux  level. 

Quickly  rotated  MRCTC 
potentiometer  to  pro- 
vide corresponding 
(♦)  1300  |y?  step  change 
in  reactivity. 


Retracted  rods  to  21,9  cm; 
allowed  flux  to  decay  to 
about  4.0  x  106. 

Rotated  MRCTC  potentiometer 
an  additional  amount  corr- 
esponding to  about  1300  pyO 
to  simulate  homogeneous  core 
conditions  and  brought  to 
steady  flux  level  of  4,5  x 

Quickly  reset  MRCTC  to  or- 
iginal setting  correspond- 
ing to  glory  hole  empty 
conditions  thus  providing 
a  step  change  in  reactiv- 
ity of  about  (-)  2600  |tx?  j 
allowed  flux  to.decay  to 
about  0,54  x  10  . 

Inserted  rods  all  the  way; 
permitted  flux  to  increase 
to  about  4,5  x  10^  and  ^ 
leveled  at  about  4.0  x  10  . 


As  indicated  in  Table  3,  the  position  to  which  simulator  rods  are  retracted 
in  each  case  is  different  from  that  given  by  the  rod  position  indicator  on 
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the  reactor  console.  This  implies  that  the  reactivity  potentiometer  con- 
tributes more  reactivity  per  centimeter  of  movement  than  do  the  actual  rods. 
By  a  trial  and  error  procedure  of  slightly  altering  the  value  of  the 
paralleling  resistance  of  the  reactivity  potentiometer,  the  actual  and 
simulator  rod  positions  might  be  more  closely  brought  into  agreement. 
The  sequence  of  events  describing  Figure  13  is  given  in  Table  4. 


TABLE  4.   SEQUENCE  OF  EVENTS  FOR  FIGURE  13 


TIME  IN 
MINUTES 


REACTOR  SEQUENCE 


CORRESPONDING  SIM- 
ULATOR  SEQUENCE 


0 

0.4 

2.1 
3o0 


Approximately  steady 
flux  level;  rods  @  15  cm; 
glory  hole  empty;  temp, 

21°C. 


Inserted  rods  at  fast 
speed. 

Rods  up;  flux  permitted 
to  increase  on  17.5 
seCo  period. 

Retracted  rods  and 
brought  to  steady  flux 

level  of  7.0  x  1060 

Scrammed  reactor  and 
permitted  flux  to  decay 
to  observe  effects  of 
delayed  neutrons 0 


Steady  flux  level; 
rods  8  10.1  cm;  MRCTC 
potentiometer  set  for 
proper  maximum  reactivity 
conditions . 

Inserted  rods  at  fast 
speed e 

Rods  up;  flux  permitted 
to  increase  on  17.5 
sec,  period. 

Retracted  rods  and  brought 
to  steady  flux  level  of 

7.0  x  106. 

Scrammed  simulator  and 
permitted  flux  to  decay 
to  observe  effects  of 
delayed  neutron  networks. 


Simulation  on  scramming  the  reactor  is  difficult  to  reproduce  because  of 
dashpot  action  on  the  rods  as  they  are  ejected  from  the  core.  This  action 
was  approximated  by  manually  rotating  the  reactivity  potentiometer  to 
the  all- rods-out  position0  As  a  result,  the  simulator  decay  curve  on 
scramming  appears  to  be  steeper  than  the  reactor  curve.  However,  the 
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general  nature  of  the  two  decay  schemes  is  quite  similar  with  the  flux 
levels  of  reactor  and  simulator  being  3.4  x  10^  and  2#3  x  io3  respectively 
four  minutes  after  scramming. 

The  sequence  of  events  describing  Figure  14  is  given  in  Table  5. 


TDfi   IN 
MINUTES 


TABUS  5.     SEQUENCE  OF  EVENTS  FOR  FIGURE  H 

REACTOR  SEQUENCE  CORRESPONDING  SIMULATOR 

SEQUENCE 


2o5 


Approximately  steady 
flux  level;  rods®  20.0 
cm;  glory  hole  empty; 
temp.  21  °C. 

Quickly  inserted  all 
poly,  graphite,  and  lead 
discs  into  glory  hole  to 
provide  a  step  change  from 
about  (-)  1300  iy?  to 
(+)  1300  jfk/o     conditions. 
Permitted  flux  to  increase 
to  about  0.73  x  106. 


Quickly  removed  all  poly 
to  provide  step  change  from 
a  supercritical  condition 
of  about  (♦)  1300  p./?     to 
a  subcritical  condition  of 
about  (-)  1300  |y0  e 


Steady  flux  level  rods 
@  22.0  cm;  MRCTC  potentio- 
meter set  for  proper  max. 
react,  conditions. 

Quickly  rotated  MRCTC  potentio- 
meter to  simulate  step  change 
of  about  2600  pp    bringing 
simulator  from  a  subcritical 
condition  of  about  (-)  1300  ty*> 
to  a  supercritical  condition 
of  about  (♦)  1300  V/O     . 
Permitted  flux  level  to  in- 
crease to  about  0.73  x  106# 

Quickly  rotated  MRCTC  potentio- 
meter back  to  original  setting 
corresponding  to  a  step  change 
in  reactivity  from  about  (+) 
1300  nz?  to  about  (-) 
1300  /o     o 


The  current  output  in  \l\l    amps,  from  an  ion  chamber  feeds  the  linear 
recorder  on  the  reactor;  the  appropriate  scales  are  labeled  in  the  three 
figures o  Corresponding  to  these  are  the  voltage  scales  on  the  simulator 
which  are  also  labeled  in  the  figures.  For  example,  the  1  x  10"   H^Ha 
scale   on  the  reactor  corresponds  to  the  100  volt  scale  on  the  simulator, 
the  1  x  10~%ua    corresponds  to  the  10  volt  scale,  the  1  x  10"10  jiya 
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FIG.  14.  FLUX-TIME  COMPARISON!    CURVE 


to  the  1  volt  scale,  and  the  1  x  10"    ^a  to  the  0.1  volt  scale.  The 
reading  in  ^a  as  well  as  the  reading  in  volts,  of  course,  are  proportional 
to  a  certain  flux  level.  The  flux  level  at  the  center  of  the  core  at  a 

power  level  of  100  raw  is  given  as  4.5  x  106,  This  corresponds  to  a 

—8 
0.45  x  10""   up.  a  output  from  the  ion  chamber  or  a  45  volt  output  from 

the  simulator;  all  other  values  are  scaled  accordingly. 

7»     Conclusions 

As  evidenced  by  the  comparison  curves  of  Figures  12  through  14,  reason- 
ably accurate  simulation  of  the  time  dependent  behavior  of  a  specific  low 
power  thermal  reactor  may  be  readily  accomplished  despite  the  many  simplifying 
assumptions  made  in  the  development  and  electrical  duplication  of  the 
governing  transient  equations.  The  final  simulator  with  its  relatively 
simple  control  system  provides  a  useful  tool  with  which  to  observe  typical 
reactor  responses  resulting  from  the  many  sequences  of  operations  to  which 
a  reactor  of  this  type  might  be  subject. 
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APPENDIX 
CIRCUIT  ANALYSIS  DETERMINATION  OF  CURRENT  Ig 


£^0      For    ANY     A 


>X 


8 


+  E 


-E 


Assume  current  directions  as  indicated.  Therefore, 


b    a      8 


Now,  let  A 
Then 


=  proportional  distance  of  P  from  mid  or  zero  position. 


-B  -  0 
R  +  P/2   (1  ♦  A) 


Therefore , 


i   = 
a 


0  -  E 
R  +  P/2   (1  -  A) 


=      -E E 

h   "  R  +  P/2   (1  ♦  A)       R  +  P/2   (1  -  A) 
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or         Ig     =    E 


=     E 


1 1 

R  ♦  P/2      (1  -  A)     "     R  ♦  P/2      (1  ♦  A) 


2  A   P/2 

£r  *  P/2      (1  -  A)] 

£R  ♦  P/2      (1  ♦  A    )} 

E  A   P 


R2   ♦  S  H  (1   ♦  A)   ♦    (1  -  AT]     ♦£!<!-  A*J 


X8     = 


A  E 


S2  *     R     ♦     |      (1  -  A*) 


or 


h  = 


A  E 


Q 


►  f)   ♦  Js   (1  -  A*) 
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APPENDIX  II 
SCALING  FACTOR  CALCULATIONS 
Referring  to  scaling  factor  equation,  namely 


(22) 


let  Rq  =  10  Nfegohms.  This  is  an  arbitrary  value  and  is  purposely  chosen 
large  to  permit  reasonable  values  of  voltage  V  which  is  dependent  on  the 
very  short  neutron  lifetime  of  this  reactor. 
Next,  assume  a  value  of  neutron  lifetime  1  of  10   sec.  Since  the  reactor 

will  not  be  far  from  critical  at  any  time,  1  and  lj  will  be  essentially 

■a- 
equal  to  each  other.  Hence,  as  a  first  approximation  10  will  also  be 

-U 
assigned  a  value  of  10   sec . 

Now  from  control  rod  calibration  data,  &        can  be  determined.  In  view 

/  max 

of  the  fact  that  normal  operating  temperatures  are  about  23°  C,/^ 
will  be  determined  at  this  temperature  assuming  a  homogeneous  core  (i.e., 
moderator,  reflector,  and  shielding  material  in  the  glory  hole).  The 
operating  range  of  reactivity  here  is 

'(-)  .034775  to  (♦)  .004550 

Therefore,  /O    can  be  taken  as  .034-775;  and  an  imposed  unbalance  in  the 
resistance  values  of  R  will  prevent  any  positive  reactivity  value  greater 
than  .004550.  In  other  words,  when  A  is  (-)  1,  /3=  " ^(mmjw) 
A  relationship  must  now  be  established  between  voltage  E  and  neutron 
population  density  n.  From  known  data, 
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0  at  100  mw  =  4.5  x  10  neuts/cm2/sec 

and  v  =  1.2?  x  10^Yt°K   cm/sec 

or   v  =  2.22  x  10   cm/sec  at  23°C. 

Therefore,  n  =  |  =  "44^-^5  =  20»25  neuts/cm3 
'      v      2.22  x  Kk 

at  rated  power  level,  100  mw. 

Because  the  maximum  output  of  the  amplifier  is  100  volts  and  because  it 


might  be  desirable  to  directly  correlate  flux  with  observed  voltage,  45 
volts  is  made  to  correspond  to  4. 5  x  10  neutroi 
constant  of  proportionality  0*  is  then  given  by 


volts  is  made  to  correspond  to  4. 5  x  10  neutrons  per  cm2  per  sec.  The 


J   =  -£-  =  -4SL  =  2  22 
w     n     20.25      # 


Next,  a  suitable  value  of  potentiometer  resistance  P  is  selected  so  that 
after  modifications  by  an  external  resistance  to  a  center  tap  connection, 
a  reasonable  value  of  total  resistance  is  still  available,  A  20K  pot 
meets  this  requirement;  and  after  modification  as  indicated  in  Figure  5t 
an  overall  resistance  of  13.33K  results   (13. OK  measured).  It  is  to  be 
noted  that  this  value  is  not  that  which  is  to  be  substituted  for  P  in 
equation  (22).  A  P  „ „  must  be  determined  which  will  be  compatible  with 
the  derivation  of  the  reactivity  potentiometer  relationships  given  in 
section  4.  As  far  as  the  circuit  is  concerned,  for  zero  or  critical  set- 
ting of  the  potentiometer  the  resistances  on  each  side  must  be  balanced. 
Therefore,  since  a  (-)  .034775  reactivity  is  required  on  one  side,  or  in 
terms  of  resistance 


V  .034775  ♦  .00455    *         \*rc?)    \±j,       xx.i  *-**• 
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This  implies  a  value  of 

Peff  =  2  x  11.5  U  =  23.0  KXL 

which  should  then  be  utilized  in  (2°)   in  place  of  P.  Now,  because  the 
actual  resistance  of  P  is  only  13.0  K,  then  10K  additional  resistance  is 
required  in  loop  2  as  indicated  in  Figure  6,  or  R2  =  R,  ♦  10K-O.  .  Values 
of  p.  are  readily  available  and  are  given  in  the  summary  of  scaling  factor 
values  at  the  end  of  section  5.  Before  calculations  begin,  a  value  of 
R  must  also  be  selected.  First,  however,  since  it  is  on  this  resistance 
that  the  individual  values  in  the  six  Re  networks  will  depend,  investiga- 
tion to  ensure  reasonable  magnitudes  for  the  Re  delay  networks  is 
necessary.  Based  upon  this  investigation,  a  final  value  of  R  =  119K  yields 
suitable  magnitudes  for  all  six  delay  components.  Therefore, 
R,  =  119K  and  R2  =  119K  +  10 K  =  129  Kjl   .  It  is  also  noted  here  that  R^ 
5P  which  yields  better  than  1^  linearity  between,©  and  A.  Calculations 
can  now  proceed  as  follows: 

R(l  ♦  B/P)/Omx   =  119  x  10*  (1  ♦  |12)  0.034775 


=  2.55  x  104jfc 


Now,  from  (22) 

JC7  J 


,*         R(l  ♦  B/P)jO 
1o  / 


max 
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In  10    * 

or     c, 


'7  R(l  *  R/PVO^  2^5  x  1Q4 


c?     =     .00392  jif 
Likewise , 


or 


J  _  _i 

r  (i  ♦  m^  '  hn± 

R  (1  *  W/°mx  _  2.55  x  10* 
Pi  Pi 


i 

Pi 

A.sec 

R.megs 

©i  »* 

lo 

.00025 

14 

102 

0.0007 

2 

.00085 

1.61 

30 

0.0207 

3 

.002a 

0.455 

10.6 

0.208 

4 

.00213 

0.154 

12.0 

0.544 

5 

.00166 

0.0315 

15.4 

2.06 

6 

.00025 

0.0125 

102 

0.785 

p     =    IE  ^     =  .00755 

Next,  the  method  for  deciding  upon  a  suitable  value  of  source  voltage, 
V,  as  discussed  in  section  5»  is  as  follows: 

B{1-K.ff)       -/°a 


30  ~~     1*  $ 
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Now,  n  ~     =  Power  ratio  x  n  -,*-> 
7     @30  jaw  @100mw 


30  x  10"6 
100  x  10~3 


x  20.25 


-3         3 
or  n  A«rt      =  6. OP  x  10   neuts  /  cm  /  sec 
@30  ^w 


Keff  *  1 
p  -    -*y •  safeties  up  position  =  (-)  .009625 

I  eff 


from  control  rod  calibration  data.  From  this  then, 


Keff      "    -99CW 


or  (1-Keff)      =    .00953 


Therefore , 

S01*  =  .579  x  10"4" 

and    S010*=  .5^5  x  10~4 


Therefore,        >,  p  ,„  *% 

v     J  y   (S°(max)  loJ 

R  (1  ♦  R/P)/2 


max 


or         v  =    (2.22) '(10  x  1Q6)'(.5?5  x  10~4) 

(2.55  x  104) 


V  =   51.0   mv 

It  is  to  be  noted  that,  in  the  case  of  the  simulation  of  the  AGN-201, 
S0  and  S0 /   \  are  the  same  since  it  is  not  desired  to  vary  the  source 
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strength  (voltage);  and  therefore  only  the  one  magnitude  determined  above 
is  needed  to  represent  the  effectively  constant  Ra-Be  neutron  source  in 
the  reactor. 

Finally,  if  1  is  taken  as  10   sec,  then  the  value  of  S0  is  readily 
found  to  be 


-4 
s  =  «^79  S  10    =   ,579  neuts/cm3/sec 

Hf4 
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